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Abstract: We examine the Non-Universal Higgs Model (NUHM) with the help of a Markov 
Chain Monte Carlo scanning technique and Bayesian statistics. The method allows us to reveal 
more fully the remarkably rich and complex structure of the model and also some new interesting 
features. In particular we find that, for large ranges of parameters, the neutralino dark matter in 
the model is higgsino-like with mass close to 1 TeV. This feature does depend to some extent on 
the prior but generically results from a mild focussing effect in the running of one of the Higgs 
masses that we also point out. On the other hand, the usual bino-like neutralino is also present, 
and in some sense is favored by a more robust prior and a better average chi-square. In spite of 
experimental constraints often favoring different regions of parameter space than in the Constrained 
MSSM, most observational consequences appear fairly similar, which will make it challenging to 
experimentally distinguish the two models. 
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1. Introduction 

Softly broken low-energy supersymmetry (SUSY) has many attractive features j^]. For 
example, unlike the Standard Model (SM), it provides an elegant solution to the gauge 
hierarchy problem and a natural weakly-interacting dark matter (DM) candidate, in addi- 
tion to accounting for gauge coupling unification. On the other hand, SUSY itself has to 
be (softly) broken in order to make contact with reality, which in the general Minimal Su- 
persymmetric Standard Model (MSSM) introduces a large number of new free parameters, 
namely the soft masses. Because of SUSY's natural link with grand unification theories 
(GUTs), one often explores SUSY models by imposing various boundary conditions at the 
GUT scale. The most popular model of this class is the Constrained MSSM (CMSSM) @, 
in which not only do gaugino soft masses unify to mi/2 but also soft masses of all the 
sfermions and Higgs doublets unify to ttt-q. These parameters, along with a common tri- 
linear mass parameter and the ratio of Higgs vacuum expectation values tan /?, form 
the four continuous parameters of the CMSSM. The relative simplicity of the model makes 
a very attractive playground for many studies. 

On the other hand, precisely because of its economy, the CMSSM may be missing some 
features of unified models with less restrictive boundary conditions at the unification scale. 
In particular, the assumption of Higgs (soft) mass unification with those of the sfermions 
does not seem strongly motivated since the Higgs and matter fields belong to different 
supermultiplets. One explicit example where this is the case is a minimal SO{10) super- 
symmetric model (MSOioSM) |^], which is well motivated and opens up a qualitatively new 
region of parameter space [Q]. Models like this provide a good motivation for exploring a 
wider class of phenomenological models in which the soft masses rrin^ and rnn^ (as defined 
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at the GUT scale) of the two Higgs doublets are treated as independent parameters and 
which come under the name of the Non-Universal Higgs Model (NUHM) [|5|. ^ 
In the NUHM, there are therefore six continuous free parameters: 

mo,mi/2,tanP,Ao,mH^ and niHa- (1-1) 

The Renormalization Group Equations (RGEs) are then used to evaluate masses and cou- 
plings at the electroweak scale and the Higgs potential is minimized in the usual way. 
Electroweak Symmetry Breaking (EWSB) conditions for NUHM read: 

„ m?r — 171% tan^ /? 1 n 
^ tan2/3-l 2 ^ ^ ^ 

"^A = mjj^ + mjj^ + 2^^ (1.3) 

where rriA stands for the mass of the pseudoscalar Higgs A and /U is the SUSY-preserving 
Higgs/higgsino mass parameter. In the above equations all the parameters are evaluated 
at the usual electroweak scale MsusY = ^m^^ m^^ (where m^^ j-^ denote the masses of the 
scalar partners of the top quark), chosen so as to minimize higher order loop corrections. 
At MsusY the (1-loop corrected) conditions of electroweak symmetry breaking (EWSB) 
are imposed and the SUSY spectrum is computed at mz- 

Like in the CMSSM the sign of fi remains undetermined. On the other hand, unlike 
in the CMSSM, because of the larger number of free parameters, both /i and uia can 
be treated as free parameters in place of high scale parameters uiHu and mn^- This 
will have an important impact on the properties of the model. In particular, the ability 
to effectively choose the position of the A funnel or to tune fi for a given point in the 
parameter space to give a correct relic abundance of the neutralino DM will lead to very 
different phenomenological predictions, as we will see below. 

The moderately increased number of free parameters of the NUHM has been shown 
to lead to a rich and distinct phenomenology (see for example |8|, ^, |lO|, |ll| and references 
therein). In particular, there is a larger variety of choices for the lightest superpartner 
(LSP). In the CMSSM the LSP is either the lightest neutralino or the lighter stau or, in some 
relatively rare case, the lighter stop. In contrast, in the NUHM, the LSP can in addition 
be a sneutrino or right handed selectron Q. Assuming the LSP to be the dark matter in 
the Universe, eliminates states that are not electrically neutral and leads to a non-trivial 
constraint on the parameter space. However, the near degeneracy of many states with the 
LSP leads to a great variety of co-annihilation channels. Also, given that tha can now be 
treated as effectively a free parameter, the resonance channel can be important in different 
ways from the CMSSM. Also, as we will see, in a large region of the parameter space the 
neutralino will actually be higgsino-like but giving the correct dark matter abundance. And 
although in general a lighter higgsino-dominated LSP underproduces dark matter as 
its mass increases the transition to higgsino dark matter can provide an acceptable relic 
density. In general these effects can also conspire to give the correct relic density in different 
parts of the parameter space. 



reduced version of NUHM with mn,, = has also been examined in several papers, eg 
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The larger number of parameters makes a full exploration of the NUHM parameter 
space even more challenging than that of the CMSSM. Additionally, varying relevant Stan- 
dard Model parameters can lead to important consequences, which has been shown in 
recent studies of the CMSSM with the help of Markov Chain Monte Carlo (MCMC) scan- 
ning technique coupled to Bayesian statistics. There are distinct advantages in doing this, 
such as being able to efficiently explore the parameter space varying all inputs simulta- 
neously and being able to incorporate errors, both theoretical and experimental, as well 
as relevant SM parameters correctly. This will allow us to efficiently explore regions of 
the model's parameters that were not easily accessible in the usual fixed-grid scans, while 
addressing, at the same time, questions about the dependence of the results on priors. 

With this in mind in this paper we apply the same approach to investigate the NUHM. 
We employ the package SuperBayeS ||l^ and include all the relevant constraints coming 
from experiment to our scan. In particular, included are the branching ratio of 6 ^ sj, 
the difference 6{g — 2)^ between the experimental and SM values of the magnetic moment 
of the muon, the LEP limits on sparticle and Higgs masses, the 5 year WMAP limits on 
the relic density r^cDM^^ and several other measured but imprecisely known quantities of 
the SM that SUSY can contribute to. A full list of constraints used and the exact numbers 
used in the analysis will be given below. 

In this paper, we present several features of the NUHM which arose in light of our 
global MCMC scans. It should be noted that we take a wider range of parameters than 
previous scans in the literature (up to 4 TeV in the case of the soft masses) , and the MCMC 
technique allows us to vary all our parameters such that we also consider Aq ^ 0, while at 
the same time varying relevant SM parameters. 

The paper is organised as follows. In sec. |2|we summarise the statistical formalism that 
we employ and list the constraints that we apply. Then, in sec. |3| we present the ranges 
of NUHM parameters that are favoured by two distinctively different prior choices and 
discuss their main features, including the possibility of higgsino DM and a mild focussing 
effect. Next, in sec. ^ we present some implications for phenomenology, including direct 
and indirect SUSY searches in colliders, while in sec. |5| we discuss prospects of direct and 
indirect detection of the neutralino dark matter in the model. We conclude and summarize 
our results in sec. ^. 

2. Outline of the statistical treatment 

In comparison to earlier analyses of the CMSSM 
parameter set, defined by 

9 = {mi/2,mo, Ao, tan (3, mH^,mHa), (2.1) 

where again we fix sgn(^)= As the relevant SM parameters, when varied over their 
experimental ranges, have impact on the observable quantities, fixing them at their central 
values would lead to inaccurate results. Instead, here we incorporate them explicitly as 
free parameters (which are then constrained using their measured values), which we call 



14, nB, nq] here we have a larger base 
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nuisance parameters tp, where 



^|J = {Mt,mb{mb)^''', aUmz)^'', a,(mz)*^^). (2.2) 

In eq. (|2.2D denotes the pole top quark mass, while the other three parameters: 
mb{mb)^^ - the bottom quark mass evaluated at nib, aem{Mz)'^^ and as{Mz)^^^ - 
respectively the electromagnetic and the strong coupling constants evaluated at the Z pole 
mass mz - are all computed in the MS scheme. Using notation consistent with previous 
analyses we define our now ten basis parameters as 

m = (6l,V) (2.3) 

which we will be scanning simultaneously with the MCMC technique. For each choice of 
m various colliders or cosmological observables will be calculated. These derived variables 
are denoted by ^ = (^1,^2, • • •)> which are then compared with the relevant available data 
d. The quantity we are interested in is the posterior probability density function (or simply 
posterior) p{m\d) which gives the probability of the parameters after the constraints coming 
from the data have been applied. The posterior follows from Bayes' theorem, 

where p{d\(^), taken as a function of ^ for fixed data d, is called the likelihood (where the 
dependence of ^(m) is understood). The likelihood is the quantity that compares the data 
with the derived observables. 7r(m) is the prior which encodes our state of knowledge of 
the parameters before comparison with the data. This state of knowledge is then updated 
by the likelihood to give us the posterior. p{d) is called the evidence or model likelihood, 
and in our analysis can be treated as a normalisation and hence is ignored subsequently. 



(See instead [jl^l for an example of how the evidence can be used for model comparison 
purposes.) 

As in previous papers, we have chosen the ranges of the base parameters as follows: 
50 GeV < mo,mi/2 < 4 TeV, \Aq\ < 7 TeV and 2 < tan/3 < 62. For the Higgs soft mass 
parameters we take < mHui^^Hd < 4 TeV. Notice that by this range of mn^ and mH^-, 
and taking sgn(/i) > 0, we automatically satisfy the GUT stability constraint of ref. 
We scan through the above parameters assuming two different choices of priors: 

• flat priors in all the NUHM parameters; 

• log priors, that are flat in logmi/21 log^zo, logm/^^ and logm//^, while for Aq and 
tan (j we keep flat priors. 

As before [^], our rationale for this choice of priors is that they are distinctively 
different. The reason why we apply different priors to soft mass parameters only is that 
they play a dominant role in the determination of the masses of the superpartners and 
Higgs bosons. As we shall show below, the choice of log priors is actually more suitable to 
the exploration of the parameter space for both physical and statistical reasons. From the 
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A/fpa n va 1 n p T In rPT"t,?i i n t,v 
;U (T (exper.) 


ref. 


Mt _ 


172.6 GeV 1.4 GeV 
4.20 GeV 0.07 GeV 
0.1176 0.002 
127.955 0.018 


1 
1 
1 
1 


23 
24 
24 
24 





Table 1: Experimental mean fi and standard deviation a adopted for the likelihood function for 
SM (nuisance) parameters, assumed to be described by a Gaussian distribution. 



physical point of view, log priors explore in much greater detail the low-mass region, which 
exhibits many fine-tuned points that can easily be missed by a flat prior scan. From the 
statistical point of view, log priors give the same a priori weight to all orders of magnitude 
in the masses, and thus appear to be less biased to giving larger statistical a priori weights 
to the large mass region, which under a flat prior has a much larger volume in parameter 
space, similar to the CMSSM case |19|. We will show below how this leads to a large 
"volume effect" in the posterior, which does not reflect the average quality of fit of those 
points. In fact, here the volume effect is expected to be even larger than in the CMSSM, 
due to the larger number of parameters. For these reasons, we regard the results obtained 
with the log prior as the more robust and we present flat prior cases mostly for comparison. 

The SM parameters are assigned flat priors and are then constrained by applying 
Gaussian likelihoods representing the experimental observations (see table |l|). The issue of 
priors is irrelevant for SM parameters, as they are directly constrained by observation. The 
predictions for the observable quantities are obtained by using SoftSusy 2.0.5 and DarkSusy 
4.0 |2^, ^] as implemented in the SuperBayeS code. The likelihoods for the relevant 
observables are taken as Gaussian (for measured observables) with mean fi, experimental 
errors a and theoretical errors r (see the detailed explanation in [^). In the case where 
there only an experimental limit is available, this is given, along with the theoretical error. 
The smearing out of bounds and combination of experimental and theoretical errors is 
handled in an identical manner to |15|. 

As ever, any points that fail to provide radiative EWSB, give us tachyonic sleptons 
or provide the LSP which is not the lightest neutralino are excluded. As in previous 
works Ip^ , 18 1, we adopt a Metropolis-Hastings MCMC algorithm to sample the pa- 
rameter space. We have also cross-checked our results by employing the more recently 
implemented MultiNest algorithm |22, 19 1 and the findings are compatible (up to numeri- 
cal noise). The results presented in the rest of this paper are obtained from the 25 chains 
that were used, garnering a total of 3 x 10^ samples each, with an acceptance rate of around 
4%. Convergence criteria are the same as in our previous papers |15, |l^, 18]. 

Although in this paper the profile likelihood is not used (see instead |19| for a detailed 
analysis), some best fit points will be shown in the 2D posterior plots for comparison. 
This is useful to identify volume effects arising from the prior, namely situations where the 
density of posterior samples (and hence, the corresponding posterior probability density) 
is large due to a large volume of parameter space and not because of large likelihood. 
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Observable 



mw 
siv? 6*2^ 

6{g -2)^ X 10^0 
BR{B Xs-i) X 10^ 

BR{Bu Tu) X 10^ 



rrih 

sparticle masses 



Mean value 



Uncertainties 
a (exper.) r (theor.) 



80.392 GeV 
0.23153 
27.5 
3.55 

17.33 ps-i 

1.32 

0.1099 



29 MeV 
16 X 10"^ 
8.4 
0.26 

0.12 ps-i 

0.49 

0.0062 



15 MeV 
15 X 10"^ 
1 

0.21 

4.8 ps-i 
0.38 

0.117^/1^ 



Limit (95% CL) 



r (theor.) 



< 5.8 X 10-s 

> 114.4 GeV (91.0 GeV) 

f{mi) 



14% 
3 GeV 
negligible 



See table 4 in ref. |15| 



ref. 



m 
m 



ref. 



11 
1311 

m 



Table 2: Summary of the observables used in the analysis. Upper part: Observables for which a 
positive measurement has been made. 5{g — 2)^ denotes the discrepancy between the experimental 
value and the SM prediction of the anomalous magnetic moment of the muon [g — 2)^. For central 
values of the SM input parameters used here, the SM value of BR{B Xg-f) is 3.11 x 10"'*, while 
the theoretical error of 0.21 x 10~^ includes uncertainties other than the parametric dependence 
on the SM nuisance parameters, especially on Mt and as{Mz)^^ ■ For each quantity we use a 
likelihood function with mean fi and standard deviation s = \/ + t^, where cr is the experimen- 



tal uncertainty and r represents our estimate of the theoretical uncertainty (see |15 for details). 
Lower part: Observables for which only limits currently exist. The likelihood function is given 
in ref. ||l^, including in particular a smearing out of experimental errors and limits to include an 
appropriate theoretical uncertainty in the observables. rn/i stands for the light Higgs mass while 
= g"^ {hZ Z)mssm/ g"^ {hZ Z)sM, where g stands for the Higgs couphng to the Z and W gauge 
boson pairs. 



3. Probability maps of NUHM parameters and observables 

In this section we present some results from our global scans of the NUHM parameter 
space. To start with, in fig. |l] we plot joint 2D probability density functions (pdfs) for 
some combinations of the NUHM base parameters in the log prior case. First, in the left 
panel in the plane spanned by mi/2 and mo we can see that the 68% total probability 
region (inner contour) is remarkably well confined to mostly a fairly low mass region of 
^1/2 ^ 1 TeV and mo ^1.4 TeV, where also the best fit points, marked by green triangles 
are located - we will discuss them below. Turning to the middle panel, we can see a 
preference for moderately large tan/3 < 40, as well as for positive ^O; although zero or 
negative values of Aq are not excluded. In some sense the latter may even be favoured as 
indicated by the best fit points. Finally, the new parameters beyond the CMSSM exhibit 
interesting behaviour and, as we shall see, a strong prior dependence. With the log prior 
rriH^ is fairly poorly constrained while favors rather low values. We will examine m//^ 
and in more detail below. 

Coming back to mx/2 and mo, we can see that, while the 68% total probability region 
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Figure 1: The 2D relative probability density functions in the planes spanned by three pairs of 
NUHM parameters: {1^1/2,1^0), (tan/3, ^0) a-nd {rnHj^jtriH^) for the log prior and for /i > 0. The 
pdfs are normalised to unity at their peak. The inner (outer) blue solid contours delimit regions 
encompassing 68% and 95% of the total probability, respectively. All other basis parameters, both 
NUHM and SM ones, in each plane have been marginalized over (i.e., integrated out). Green 
triangles denote the three best fitting points. 




Relative probability density 



Figure 2: The same as in fig. |lj but for the flat prior. As explained in the text, the considerable 
shift in the probability density is largely due to a volume effect from the flat prior arising from the 
large number of samples at large mass for this choice of priors. The results obtained with log prior 
(fig. ^) are to be considered more robust as they are the ones where the probability density is in 
better agreement with the average quality of fit. 



is well constrained, the 95% region (outer contour) is much wider and extends to much 
larger ranges of both parameters. This signals that the constraining power of the data is 
not sufficient to strongly confine the posterior. In order to examine this, in fig. ^ we plot 
the same quantities as in fig. || but assuming the flat prior. We find that the region of large 
1^1/2 ^ 2 TeV, as well as of 0.5 TeV < mo < 3 TeV (with a preference for mi/2 > mo) 
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is in this case much more prominent and actually favored from the point of view of the 
posterior. As we shall show below, this is actually largely due to a volume effect in the 
prior, for in fact the average the large mass region is relatively poorer than in the 

small mass region. The upper bound on the high probability regions of ?ni/2 and niQ mostly 
results from the adopted upper limit on mn^^ ■ The best fit points have also shifted upwards 
but not nearly as much and remain far away from the 68% total probability region. This 
suggests that the resulting picture suffers from the previously mentioned volume eff'ect due 
to integrating over the additional dimensions: at large the fit to the data is moderate 
but there is a large volume of such points.^ In Bayesian statistics such situations appear to 
be be favored because of the large number of samples contained in the large volume implied 
by the choice of measure contained in the prior. Thus the flat prior appears to override 
the preference for the small mass region (where the average is generally better), in the 
sense that it imposes a measure on parameter space that does not adequately explore the 
low mass region. For this reason below we will present results for that choice mostly in 
order to illustrate prior dependence, rather than by treating it on equal footing with the log 
prior. It should be kept in mind that the log prior results are the ones where the posterior 
distribution more fairly represents the actual quality of flt of the likelihood. It is expected 
that with better, more constraining data this prior dependence will be resolved and the 
posterior will become largely prior independent. 

This approach is further justified by comparing the other two combinations of base 
parameters presented in figs. |^ and ^. In particular, with the fiat prior the range of 
has much expanded, especially towards larger positive values. The 68% and 95% total 
probability ranges of tan /? have actually shrunk somewhat but the best fit points have 
moved from around 20 in the log prior case to nearly 40. The favored ranges of the Higgs soft 
masses show an even stronger volatility with prior choice. While, as we have noted above, 
for the log prior tuh^ is rather strongly pushed toward lower values, well below 1 TeV, 
but mn^ remains fairly unconstrained. In contrast, with the flat prior mn^ shows strong 
preference for large values close to its upper limit while rriH^ is generally unconstrained. A 
proper understanding of this seemingly uncontrollable behavior will require some further 
investigation that we will conduct below. Nevertheless, this point illustrates the difficulty 
in reliably constraining both mij„ and rnn^ with currently available data. On the other 
hand, the best fit points for both priors tend to agree more with the log prior posterior, 
which is another sign that this choice of prior leads to a posterior being better aligned with 
the high likelihood regions. 

It is worth comparing the high probability regions of the NUHM with those for the 
analogous parameters in the CMSSM. The left two panels in figs. |l] and || are directly 
comparable with the corresponding panels in fig. 13 of |jl^. The flrst point to remember 

^To estimate the magnitude of this volume effect, recaU that, for the flat prior, the volume encompassed 
by the mass range between e.g. 1 and 4 TeV is a factor 10* larger than between 100 and 400 GeV. For 
comparison, under the log prior the ratio of the volumes of the two regions is unity. Therefore, in order to 
completely override the prior volume via the likelihood, one would need at least a ~ 4.3(t preference for the 
low mass region. Clearly, current data is not constraining enough to override the prior preference for large 
mass in this case. 
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is that prior dependence was already rather strong in the CMSSM and it should come as 
no surprise that in the NUHM, with two more free parameters, the 2D joint pdfs with all 
the other parameters marginalized over, show an even stronger dependence. Interestingly, 
in both models different regions of the corresponding parameters are typically favored. In 
the log prior case, in both the CMSSM and the NUHM, rather low "^-o ranges 

are favored (a reflection of log prior's preference for low masses), but actually different 
ones and selected by different physical mechanisms. In the CMSSM this is mostly the 
the neutralino-stau coannihilation region of mo <C Tni/2 ^ 1 TeV [|l9| (plus a tiny vertical 
region of Z and h pole annihilation), while in the NUHM the analogous (mentioned above) 
"low-mass" ranges mi/2 ^ 1 TeV and mo ^1-4 TeV corresponds to the "bulk region" 
where A funnel annihilation play a dominant role. In both models, the favored ranges 
(and also best fit points) of tan (3 are rather moderate (although in the CMSSM very large 
values of around 55 are not excluded). Finally, ^40 in both models shows a mild preference 
for positive values. 

In contrast, for the flat prior the difference between the overlapping CMSSM and 
NUHM parameters is much larger, which, again, is a reflection of the volume effect in 
both cases. In the CMSSM the focus point region (FP) ||3^ of large mo ^"^1/2 is favored 
(although with the stau coannihilation region still remaining allowed) fl^, |ig|]. In contrast, 
in the NUHM a very different region of large mi/2 ^ w-o is favored, as discussed above. 
Apart from the statistical reasons given above, the underlying physics is also different. 
In the CMSSM, the FP region is to a large extent favored by an interplay between the 
current SM prediction for and an experimental world average of BR[B Xg^) |ll8| , |l9|| . 
The case of the NUHM will be discussed shortly. Also, the strong preference for very large 
tan /3 ~ 50 — 55 in the CMSSM should be contrasted with significantly lower values of the 
parameter being favored in the NUHM. This may be because in the NUHM the large tan f3 
regime affects strongly the running of the slepton masses via the Yukawa couplings which 
appear in the RG equations, making 1^,6 large (and making Yt small) which could lead to 
problems with EWSB or tachyonic sleptons | |10| |. 

In the NUHM with the flat prior, the high probability region of large mi/2 iii fig- @ 
corresponds to the (lightest) neutralino x being mostly higgsino. This can be seen in the 
right panel of fig. ^ where we plot the values of the x gaugino fraction Zg = Zf^ + Z\2 
in the m^/2)"^o plane, where Z^i is the bino fraction of the neutralino and wino 
fraction. We show three distinct regions: red dots correspond to a mostly higgsino state, 
Zg < 0.3, green squares to a mixed state (0.3 < Zg < 0.7) and blue diamonds to mostly 
gaugino neutralino, Zg > 0.7. This is a new feature of the NUHM which distinguishes 
it from the CMSSM where the neutralino is mostly a bino (also in the FP where the 
higgsino component, while larger, still remains subdominant).^ Although in general a 
lighter higgsino-dominated LSP underproduces dark matter, this is not the case for heavy 
enough higgsinos. The significance of the higgsino LSP region clearly depends strongly on 
the assumed prior, which is evident by comparing the left and right panels of fig. ^ (and the 
corresponding left panels of figs. |l] and |2|), but its presence is an interesting new feature of 

^The possibility of higgsino-Iike LSP in the NUHM has also been noticed in but not explored in 
more detail. 
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Figure 3: Values of the gaugino fraction Zg = Zf^ + ZI2 in the plane of (?7ii/2i"^o) for samples 
uniformly selected from our MC chains obtained assuming the log prior (left panel) and the flat 
prior (right panel). The color coding is as follows: red dots correspond to Zg < 0.3 (mostly 
higgsino), green squares to 0.3 < Zg < 0.7 and blue diamonds to Zg > 0.7 (mostly gaugino). The 
triangles denote the best fit point for each cloud of samples of a given respective gaugino fraction 
(of corresponding color) taken separately. The overall best-fit is for both priors in the gaugino-like 
DM region. 




Figure 4: For the neutralino mass and its gaugino fraction Zg we show the posterior for the 
log prior (left panel) and flat prior (right panel). 



the phenomenology of the NUHM that is not present in the CMSSM. More quantitatively, 
we find for the log prior that the probability of the higgsino region is ~ 12%, which increases 



-10- 



to ~ 70% for the flat prior case. Although the latter figure is strongly influenced by the 
volume effect, even in the more conservative case of a log prior we can still conclude that 
the higgsino DM region has a sizable probability in the NUHM (of order 10%). We can also 
investigate the beahviour of the for each of the three regions separately and for both 
priors. We find that the overall best-fit is always to be found in the region where the LSP 
is mostly gaugino (blue samples in fig. for both priors. However, the average in each 
of the three clouds of samples is almost constant, which clearly shows that volume effects 
are at work and that the relative population of samples in the three regions is strongly 
influenced by the choice of prior. ^ 

In order to examine in more detail the higgsino DM in the NUHM, in fig. ^ we present 
a 2D joint pdf in the neutralino mass and its gaugino fraction. The higgsino region is 
clearly pronounced at m^^ ~ 1 TeV in the flat prior case but it is much less pronounced 
in the log prior case where the neutralino is mostly a bino. This will affect the ensuing 
predictions for DM searches, as we shall see later. 
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Figure 5: The running of m*^^{Q) (left panel) and ra*u^{Q) (right panel) with the log of energy 
scale t = In Q. We take = 4 TeV, mi/2 = 3 TeV, tan/3 = 5 and, in the NUHM case (red sohd 
lines), mo = 2 TeV. In order to facilitate comparison with the CMSSM (blue dash-dotted lines), for 
each curve we initially set — "^h^ and then evolve differently for each model. In the CMSSM, 
for each curve we take toq — rriu^ — itih^ , which is why the middle curves overlap. 

For further discussion it will be convenient to introduce the running parameters rn^^ (Q) 
and 'rn*jj^{Q), where Q is the energy scale, defined as 

^h.JQ) = K„,.(Q)) \<.JQ)\'^'- (3-1) 

*We have also investigated the impact of the constraint coming from the anomalous magnetic moment 
of the muon, 5{g — 2)^, which has been shown to be in some tension with fe — > 37 [^. If one removes 
5{g — 2)p from the analysis, the average ^nd the best-fit values for the three different region become very 
close to each other. This shows that the fit is being driven to a large extent by the (somewhat controversial) 
5{g — 2)p constraint. As such, one has to exercise some care in interpreting the ensuing conclusions. 
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Since in the RGEs the running Higgs soft mass parameters (Q) and rriH^ (Q) appear 
only in squares, which can become negative, the parameters mj^ ^^^^ convenient to 
deal with in the sense that they adequately reflect both the magnitude and the sign of 
the respective parameters m'j^ d^^^' particular, the quantities ^, without any 
arguments, will denote the respective running quantities evaluated at Q = MsusY; 

"^?^.,. = "^k,.(Q = ^susY). (3.2) 

We are now in the position to address the apparent strong prior dependence of 
and ruHd in figs. |l| and This feature, and also the presence of the higgsino LSP, originate 
from another new interesting feature of the NUHM parameter space which is the existence 
of a mild focussing effect in the RG running of m//^, akin to that in the CMSSM |^ . 
The effect is illustrated in fig. ^ where red solid lines show the running of the NUHM 
parameters m*jj^{Q) (left panel) and rn*j^^{Q) (right panel) with the log of energy scale 
t = In Q. For each case we take Aq = 4 TeV, mi/2 = 3 TeV and tan/3 = 5. In order to 
facilitate comparison with the CMSSM (blue dash-dotted lines), for each curve we initially 
set TnH^ = nT'Hu and then evolve differently in each model. In the NUHM we fix also 
mo = 2 TeV while in the CMSSM, for each curve we take ttt-q = TriHa — ^Hui which is 
why the curves in the middle case overlap. One can see that, in the NUHM the running of 
is stronger for larger GUT values of the parameter, but it is not as strong as in the 
CMSSM. 

As a result, we can see some "squeezing" of m,^^ compared to the GUT values rriHy,, 
while this is not the case with the soft mass parameter. This is shown in fig. |6| where, in 
the flat prior case (right panel) the fairly narrow inner 68% total probability region of rn'^^ 
between some —1 TeV and —0.5 TeV corresponds to large values of rnn^ (compare fig. |2|), 
close to the assumed upper limit of the prior, nearly independently of the Ha soft mass 
parameter. It is clear from this that ni'^^ is to a large extent constrained by the focussing 
effect in the RG running, while the most probable ranges of n^n^ si's more prior dependent. 
On the other hand, in the log prior case prefers to be smaller than for the other prior 
(compare right panels of figs. || and ||) basically in order to arrive at bino-like neutralino 
with correct relic density, as we will show below. In this case, however, mn^ is confined to 
preferably fairly small values (< 2 TeV) as otherwise one ends up with tachyonic sleptons. 
(We can also again see that, in the log prior case the high probability region agrees rather 
well with the best fit points, while this is not so with the flat priors.) The corollary to this 
is that, in general we have more freedom in obtaining a phenomenologically desired range 
of values at the electroweak scale by appropriately choosing ruHy, at the GUT scale. Clearly 
that is not so in the CMSSM as one can never attain smaller fi here due to the stronger 
FP behaviour essentially focusing to nearly one point at the EW scale. The two distinct 
branches visible in fig. |6| correspond to two distinct neutralino regimes, as shown in fig. |^. 
In the horizontal branch the neutralino is mostly a higgsino (red dots) while the other, an 
inverted 'V shaped region around = gives us a mostly bino (blue diamonds). 

Above the reach of current collider limits, the high probability regions of mn^ , nT-H^ and 
the other soft parameters are primarily determined by requiring the correct dark matter 
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Figure 6: For to"^^ and fn""^^ we show the posterior for the log prior (left panel) and fiat prior 
(right panel). 
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Figure 7: Values of the gaugino fraction Zg in the plane of {m'^^ , m*^^ ) for samples uniformly 
selected from our MC chains obtained assuming the log prior (left panel) and the flat prior (right 
panel). The color coding is as follows: red dots correspond to Zg < 0.3 (mostly higgsino), green 
squares to 0.3 < Zg < 0.7 and blue diamonds to Zg > 0.7 (mostly gaugino). 



abundance.^ It is well known that, in order to satisfy this constraint, the LSP neutralino 
can either be mostly bino-like (hke in the CMSSM) if the bino soft mass Mi < or 
else a sufficiently heavy higgsino-like state with < Mi. With flat priors imposed on 
the soft masses it is the second choice that is primarily realized in the NUHM, because, 

^Exploratory runs with the constraint switched off yield much wider ranges of parameters. 
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as explained above, by starting with large enough mj^„ at the GUT scale one arrives, via 
the mild focusing effect (see fig. ^ at less negative values of at the EW scale. In 
the limit Im^J ^ ^^'^ large enough tan/3, eq. ( |1.2| ) would imply n ~ 

reality, in the NUHM this limit is often violated and as a result /x comes out somewhat 
larger. In any case, with flat priors one obtains small enough values of /x and higgsino dark 
matter. On the other hand, ~ has to be large enough to give an acceptable relic 
density, as mentioned before. Numerically, this leads to > 0.8 TeV, which translates (via 
0.4mi/2 ~ Ml > |/i|) to mi/2 > 2 TeV; see the left panel of fig. |. The second EWSB 
condition, eq. (^), then implies the approximate relation m\ ~ m|^^ + fj? which we have 
checked numerically to hold. 

The condition fi > 0.8 TeV is reflected as the left side of an inner (68%) vertical contour 
in the right panel of fig. ^ where we show the joint 2D posterior in the plane (^, m^). On 
the other side, as increases, so does fi^/i^ which quickly becomes unacceptably large. As 
a result one finds a strong concentration at /i ~ 1 TeV, corresponding to higgsino LSP. 




Figure 8: For /i and itia we show the joint 2D posterior for the log prior (left panel) and flat 
prior (right panel). 

The higgsino LSP at ~ |^| ~ 1 TeV is further emphasized in fig. ^ where we show 
values of the gaugino fraction Zg in the plane of (m-^,m^) for samples uniformly selected 
from our MC chains obtained assuming the log prior (left panel) and the flat prior (right 
panel). The color coding is as before: red dots correspond to Zg < 0.3 (mostly higgsino), 
green squares to 0.3 < Zg < 0.7 and blue diamonds to Zg > 0.7 (mostly gaugino). The 
vertical higgsino region at ~ 1 TeV is strongly pronounced for the flat prior although 
is also visible in the log prior case whose impact in pushing masses to lower values is again 
visible. 

On the other hand, the diagonal branch in fig. ^ corresponds to the second way of 
arriving at the correct relic density, namely via bino-like neutralino which is clearly also 
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Figure 9: Values of the gaugino fraction Zg in the plane of {ra-^^rriA) for samples uniformly 
selected from our MC chains obtained assuming the log prior (left panel) and the flat prior (right 
panel). The color coding is as follows: red dots correspond to Zg < 0.3 (mostly higgsino), green 
squares to 0.3 < Zg < 0.7 and blue diamonds to Zg > 0.7 (mostly gaugino). 



realized in the NUHM. In this case the LSP mass extends over a range of values and 
the correct relic density is achieved via a A resonance when ~ ^m^^. This feature can 
be seen in both panels of fig. ^. It also explains a kink at /i ~ 1 TeV in the right panel of 

In the case of the log prior, with its built-in emphasis on low masses, the preference 
for the bino and the higgsino DM regions in some sense gets interchanged, as can be seen 
by comparing the left and right panels of the above figures. The effect of the log prior is 
twofold. Firstly, it serves to put the emphasis to low mi/2 where, even if one can achieve 
fi < Ml, the higgsino will be too light to give the correct relic density. Coupled with this, 
the emphasis on the low mass region of also plays a role: if one starts with too low 
values of mn^ at the GUT scale, one arrives at too large values of |//| and an incorrect 
relic abundance. The only way to satisfy the cosmological constraint is then to rely on the 
bino-like DM which requires small enough Mi and therefore mi/2- This is precisely what 
one can see in fig. |^ and subsequent figures showing the log prior case. 

In summary, the parameters of the NUHM exhibit a rather complex structure. Clearly 
the two priors adopted here give us sizeably different results, but even the strong preference 
for the low mass region encoded in the log prior is not sufficient to complete kill off the 
much bigger regions of large mi/2 and mo which are strongly favored by the flat prior. In 
some senses we have two distinct regimes here. One has a 1 TeV higgsino neutralino and 
corresponds mostly to large values of the soft mass parameters, whilst the other is more 
similar to the situation seen in the CMSSM, with bino dark matter and a lighter spectrum. 
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4. Collider signatures 
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Figure 10: The ID relative probability densities for fl^h^ (upper left panel), BR{B Xsj) 
(upper right panel), S{g — 2)^ (lower left panel) and the light Higgs mass rrih (lower right panel). 
In each panel we show the posterior for the flat prior (dotted blue) and the log prior (long-dashed 
red) and the likelihood function (solid black). 

Now we turn to discussing some of the observables that play a key role in determining 
the favored regions of NUHM parameter space and also discuss some implications for the 
light Higgs, neutralino, chargino and gluino mass spectra. 



To start with, shown in fig. IC are the ID posterior pdfs for the four key constraints on 
the model: the neutralino relic abundance Q.-^^h'^, BR{B ^sT), the SUSY contribution 
to {g — 2)^ and the light Higgs mass m/j. It is clear that a tight WMAP constraint on fi^/i^ 
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Table 3: A table showing the values of various parameters for the best fitting point in both the 
log and flat prior case, as well as the effect of removing the {g — 2)^ constraint on the best fitting 
(i.e., log prior) point. 



in the likelihood is reasonably well matched by the shape of the ID pdfs for both priors. 
Likewise, the ID pdfs of BR{B Xg'j) for both priors show a rather strong preference 
for basically the SM value of the observable, although, naturally, with the log prior we find 
a larger tail of bigger values, much closer to the experimental average. Also, the SUSY 
contribution to — comes out rather tiny, resulting in 5{g — 2)^ ~ for both priors that 
we examined. This is an interesting characteristic feature which we have already seen in the 
case of the CMSSM |18, 19 1. Finally, the posterior for rrih shows a somewhat stronger prior 
dependence, although the case of the log prior does favor lower mh values and allows for a 
wider spread due to its built-in preference for lower values of basis parameters. Interesting, 
unlike in the CMSSM, the lightest Higgs is not necessarily SM-like and therefore the 95% 
LEP lower bound on SM Higgs mass should only be considered as indicative. We will 
investigate the Higgs sector of the model in more detail elsewhere. 

Generally, it is interesting to note that posterior pdfs for these observables do not vary 
wildly with differing prior, although the log prior in general seems to lead to somewhat 
less strongly peaked posteriors. Finally, it is remarkable that even in this rather different 



model, the pdfs for both set of priors are fairly similar to those in the CMSSM [18, |l 



In table we present best fit values for the NUHM base parameters and for a number of 
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quantities of particular interest, as well as the values. Clearly, the log prior does provide 
a better than the flat prior and in this sense the log prior does look more appropriate to 
explore the properties of the model. The best fitting point in the flat prior case is included 
for contrast and to illustrate the larger values for the base mass parameters and derived 
spectrum. Back to the log prior case, it is interesting to note that although removing the 
{g — 2)^ constraint does not spectacularly improve the fit, it does move the best fit for 
the base parameters and corresponding masses to much higher values. On the other hand, 
removing the {g — 2)^ constraint brings the best fit value for the fiat prior back in line 
with the value obtained for the log prior, thus giving a marked improvement over the case 
where it is included. That said, the main "indirect" observables Q^h?, BR{B Xg^) 
and BR{Bs (although not (^f — 2)^) remain fairly similar, in line with what is 

observed for the pdfs. 

On the other hand, the impact of the priors on superpartner masses is quite strong, 
with a whole slew of heavy particles in the fiat prior case being strongly contrasted with 
the low mass regions (and hence lighter spectra) preferred by the log prior. This point 
is illustrated in fig. |lT| where we present the ID posterior pdfs for the lightest neutralino 
(left panel) , the lighter chargino (middle panel) and the gluino (right panel) . This provides 
further evidence that the data is not yet good enough to properly constrain the model (as 
was already the case in CMSSM [|l^) and as such care must be taken with the results. In 
particular, LHC prospects for detecting the gluino, up to some 2.7 TeV, can be considered 
either very good or very poor depending on the prior. 
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Figure 11: The ID relative probability densities for the mass of the lightest neutralino (left 
panel), the lightest chargino fn^± (middle panel) and the gluino rrig (right panel). In each panel 
we show the posterior for the flat prior (dotted blue) and the log prior (long-dashed red). 



5. Dark matter signatures 

In this section we will examine implications for the detection of the lightest neutralino 
assumed to be the DM in the Universe. We will consider both direct detection via its 
elastic scatterings with targets in underground detectors, as well as indirect signatures 
of neutralino pair annihilation resulting in an additional component of diffuse gamma 
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radiation from the Galactic center and of positron flux from the Galactic halo (see [34| 
instead for a recent analysis of the detection prospects for the gamma~ray flux from nearby 
dwarf galaxies in the context of the CMSSM). The underlying formalism for both direct 
and indirect search modes can be found in several sources. (See, e.g., |36, 35 



In 



our analysis we have followed the procedure as well as hadronic matrix elements inputs as 
presented in our earlier work |15, 18, |3^. Some investigations into this area in the case of 
non-universality have also been done in the literature, see for example ref. ||^, [40| . 



To start with, in fig. 12 we present 2D posterior pdfs in the usual plane spanned by 

The left (right) panel 



the spin-independent cross section c7p^ and the neutralino mass 
corresponds to the log (flat) prior. For comparison, some of the most stringent 90% CL 



experimental upper limits are also marked |41, 42, 42, 44, 46], although they have not 



been imposed in the likelihood, as before in our studies of the CMSSM. As previously with 
superpartner masses, we can see a strong prior dependence, with the log prior generally 
favoring lower neutralino masses and therefore larger cross-sections. For the flat prior, in 
contrast, there is a well pronounced region of high probability in the higgsino LSP region of 
niy. ~ 1 TeV, as we discussed above, which is present but less emphasized in the log prior 
case. The prior dependence is certainly somewhat discomforting, although it is merely a 
reflection of the model being unconstrained by the currently available data. Nevertheless, 
the larger 95% total probability regions (outer contours) are actually more similar. In 
any case, we note that in both cases the 68% regions of total probability (inner contours) 
correspond to 0"^^ > 10"^'' pb which means that, independently of the prior, the direct DM 
detection predictions of the model will be basically completely explored by current and 
especially future 1-tonne detectors whose sensitivity reach is likely to be of that order or 
better. 




Figure 12: For the dark matter spin-independent cross section vs. the neutralino mass 
we show the 2D relative probability density for the log prior (left panel) and the flat prior (right 
panel) . 
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Next we discuss some indirect detection modes of much current interest. We follow 
the formalism and procedure outlined in ref. |^9| which we briefly summarize here for 
completeness. In our numerical analysis we rely on DarkSusy [^] to compute the fluxes. 
To start with, we compute the total 7-ray flux produced by dark matter annihilations 
in the Galactic halo, 



<l>^{An)= dE^d<^j/dE^{E^,An), (5.1) 



where the cone Ail. is centered on the direction ip and the integration goes over the range 
of photon energies from an energy threshold £^th up to m^. The differential diffuse 7-ray 
flux arriving from a direction at an angle ip from the Galactic center (GC) is given by 



*'(^..*)-E£i^I/'pJ«'.^o). (5.2) 



where aiV is a product of the WIMP pair-annihilation cross section into a final state i times 
the pair's relative velocity and dN^^/dE^ is the differential 7-ray spectrum (including a 
branching ratio into photons) following from the state i. Here we consider contributions 
from the continuum (as opposed to photon lines coming from one loop direct neutralino 
annihilation into 77 and ^Z), resulting from cascade decays of all kinematically allowed 
final state SM fermions and combinations of gauge and Higgs bosons. The integral is 
taken along the line of sight (l.o.s.) from the detector. It is convenient to separate factors 
depending on particle physics and on halo properties by introducing the dimensionless 
quantity J(V') = (1/8.5 kpc) (0.3 GeV/cm^f J^ ^ ^ dl p^^r {I, ip)) @. The flux is further 
averaged over the solid angle AO, representing the acceptance angle of the detector, and 
one defines the quantity J{AQ) = (l/An)J^^J(ip)di}. Since we are interested in the 
Galactic center, we set ip = 0. 

The flux from the GC critically depends on the dark matter halo profile at small 
Galactic radius r where dark matter density is thought to be largest. In this analysis we 
consider a generic NFW model p3] , as done previously in pSl] , but in addition also include 



a model of Klypin, et al., in |49|. It is based on the NFW model but it is fitted to the 
data from the Milky Way and takes into account the effect of angular momentum exchange 
between baryons and dark matter, and in this sense may be considered as more applicable 
to our Galaxy. The inner radius density profile for the Klypin, et al., model is ~ /--i-S 
which is quite a bit steeper than the NFW one with r~^. Close to solar radius both models 
become quite similar. 



In fig. 13 we present our predictions for the diffuse 7-ray flux produced by dark 



matter annihilation in the Galactic center vs. the neutralino mass for the two halo 
models. We also assume a conservative energy threshold £'th = 10 GeV and AQ = 10~^sr 
to match Fermi's resolution. We show the 2D relative probability density for the log prior 
(left panel) and the flat prior (right panel). For the Klypin, et al., model both the 68% and 
the 95% total probability ranges (inner and outer contours, respectively) are shown, while 
for the NFW only the 68% regions are plotted for comparison. Fermi's reach with one 
year of data is also indicated (horizontal black/dashed line |5^]). For each prior the fluxes 
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Figure 13: The diffuse 7-ray flux <i>^ produced by dark; matter annihilation in the Galactic center 
vs. the neutralino mass for some popular halo models. We show the 2D relative probability 
density for the log prior (left panel) and the flat prior (right panel). Klypin model (right panel). 



are somewhat different, although in the same ballpark. The emerging picture is actually 
fairly similar to the CMSSM [^]. In particular, it is clear that the largest uncertainty in 
assessing Fermi's prospects for detecting a signal in this model lies in the cuspiness of the 
dark matter halo profile at small radii. If the profile is a steep as in the Klypin, et al., 
model, then Fermi's prospects in this mode are very good, unlike for less cuspy models, 
like the NFW one. 

Finally we present the NUHM's predictions for positron flux from neutralino dark 
matter annihilation in the local halo. Once produced, positrons propagate through the 
Galactic medium and their spectrum is distorted due to synchrotron radiation and inverse 
Compton scattering at high energies, bremsstrahlung and ionization at lower energies. The 
effects of positron propagation are computed following a standard procedure described 



in 1 51, 52 1, by solving numerically the diffusion- loss equation for the number density of 
positrons per unit energy drig+fde. The diffusion coefficient is parameterized as K{£) = 
Ko{3°' + e°), with Kq = 2.1 x lO^^ cm^sec^^, a = 0.6 and e = E^+/l GeV, mimicking 
re-acceleration effects. The energy loss rate is given by 6(e) = tee'^, with te = 10~-'^^sec~-'^, 
and we describe the diffusion zone (i.e., the Galaxy) as an infinite slab of height L = 4kpc, 
with free escape boundary conditions. Changes in the above positron propagation model, 
especially K(e) (see, e.g., ]53|, ^), can potentially lead to variations by a factor of 5 to 10 in 



the spectral shape at low positron energy, Eg+ < 20 GeV [54|. Most high-energy positrons, 
on the other hand, originate from the local neighborhood the size of a few kpcs p2| , 55|, 
and their flux is less dependent of the halo and propagation dynamics. In order to reduce 
the impact of solar winds and magnetosphere effects on the model's predictions, it is 
useful to consider the positron fraction, defined as $g+/($g+ + $£-), where <I>e+ is the 
positron differential flux from WIMP annihilation, while ^>e- is the background electron 
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Figure 14: Predicted positron flux fraction produced by dark matter annihilation in the Galactic 
halo vs. the positron energy E^+, assuming log priors and the boost factor BF=1, for the Klypin, 
et al., model (left panel) and the NFW model (right panel). 



flux. For background e and fluxes we follow the parametrization adopted in ref. |£2| 
from ref. |53]. 



In fig. 14 we present the positron flux fraction produced by dark matter annihilation 
in the Galactic halo vs. the positron energy for the two halo models considered above 
and for the log prior case. (The flat prior case is not qualitatively different.) Also included 
are the relevant experimental data |5G, 57 1, including the recent Pamela result. It is clear 
that supersymmetric dark matter in the models like the NUHM (and also the CMSSM pS|]) 
falls far short of reproducing the Pamela result. This would remain true even if one would 
be prepared to consider an unlikely existence of very dense local DM clumps for which 
the boost factor (BF) would be unrealistically high, ~ 10^. Although a more refined 
analysis fitting signal and background simultaneously would be required to draw more 
quantitative conclusions, it is clear that the spectrum predicted by the NUHM appears to 
have a very different energy dependence from the flux observed by Pamela. This is not 
necessarily a problem for the NUHM, and other unified SUSY models like the CMSSM, as 
long as their signal remains below the observed flux, since a more conventional astrophysical 
explanation in terms of pulsar radiation may be entirely sufficient to account for Pamela 
observations [p^]. 

6. Conclusions 

The MCMC scan and Bayesian analysis of the Non-Universal Higgs Model performed in 
this paper reveal a remarkably rich and complex structure of its parameter space. While 
the properties of the model are in some aspects fairly similar to the CMSSM, we have 
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found several interesting differences. Perhaps the biggest one is the existence of higgsino- 
hke dark matter with a mass close to 1 TeV, with a minimum probability of about ~ 10%. 
This feature seems quite robust, as our exploratory scans with wider ranges of soft masses 
have confirmed, although the details are quite strongly prior dependent. The higgsino 
dark matter results from having more parameters than in the CMSSM but also from the 
focusing effect being less strong than in that model. 

In terms of observational consequences at colliders, the NUHM actually appears rather 
similar to the CMSSM, which will make it difficult to experimentally distinguish the two 
models. Again, the best prospects may be provided by finding in direct detection searches 



a ~ 1 TeV dark matter WIMP, since such a case in the CMSSM is highly unlikely [18, 19 1. 
Fermi's prospects for probing the model strongly depend on the cuspiness of the dark 
matter profile in the Galactic center while positrons produced in dark matter annihilation 
remain well below the Pamela result. 
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